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INTRODUCTION 


The core of a nuclear reactor constitutes an intense source of gamma radiation 
which gives rise to serious shielding and heating problems. To solve these problems 
one needs reliable methods for calculating the gamma ray flux and dose rate out- 
side the reactor core. 

It is, however, not generally possible to solve gamma ray problems exactly, 
although the differential equations can be formulated rather easily and the physical 
constants are fairly well known. Even simple problems tend to overtax the capa- 
bilities of the largest computing machines and practical problems can only be 
dealt with approximatively after strong simplifications. 

In this work we investigate the effects of various approximations which permit 
the calculation of gamma radiation from a reactor core by manual methods. We 
consider a heterogeneous thermal reactor with cylindrical rod fuel elements forming 
a cylindrical core. A reactor of this type is the early Swedish R3b project, which 
forms the basis for our calculations. 


We eventually conclude, that one of the methods is best suitable for computing 
machine programming. Such a program has in fact been described in a report by 
the author (1959 b). 


1. BASIC DATA OF THE R3b PROJECT 


The R3b is a DgO-cooled and DgO-moderated natural uranium reactor, designed 
to be operated at 94 MW power. The 112 fuel elements are arranged on a square 
lattice with a 25 cm pitch so that they form an approximately cylindrical core. 
Each fuel element consists of 27 cylindrical Zrecanned UOg rods enclosed in a Zr 
pressure tube. The material proportions of the core and the main dimensions are 
listed in table 1. 


Table 1. 


Basic data of the R3b project. 


Diameter of the core (2R,) 319 cm 
Extrapolated diameter of the core (2R) 388.6 cm 
Height of the core (Hj) 332 cm 
Extrapolated height of the core (H) 379.3 cm 
Diameter of fuel element 11.16 cm 
Volume of UO, in fuel element 1.118-104 cm? 
Mass " UO, " core 12.77 ¢t 
° 26.76 
yet 4.65 t 
Average density of fuel element 5.230 g/cm? 
Average power density in the core (Py) 3.543 W/em® 


The gamma sources and energy spectra of reactors of the type considered have | 
recently been discussed in a publication by the author (1959 a). Since we in this 
work, however, are dealing with manual methods, requiring the energy spectra to 
be divided into a few rough groups, it is sufficient to use the old spectra of Braun 
(1957) which also were used in the original report by Braun and Roos (1957). 
Table 2 lists the data of Braun's spectra necessary in this work. Here N,o(E) is 
the total energy release in the energy group of mean energy E, and Nesc(B) is 
the part of the total energy which escapes from a fuel element. The fraction 
NesclE)/Nioy(E) is the escape probability, a function of ru.,, where r is the radius 
of the fuel element and pe, is the energy absorption coefficient. The escape prob- 
ability has been tabulated by Placzek & al. (1953). [ 


Table 2. 


Energy groups of the spectra of Braun. 


Group interval | Mean energy E Nrot (2) Nese (B) 
(MeV) (MeV) (Mev (Mev 
O-- 3 2 17.56 6.08 
5 4 5.95 2.62 
5 = 8 7 1.59 0.62 
0-4 3 7.62 
4- 8 6 1.70 
O-- 8 25.10 9.32 


x) We use f as short-hand notation for “fissions". 


In practical shielding calculations the gamma ray sources have to be expressed 
in units of MeV/s rather than in MeVA, which is convenient for spectral re- 
presentations. Depending on the geometrical shape of the source one defines a point 
source strength S(MeV/), a line source strength SL (MeV/cm s), a surface or plane 
source strength S A (MeV /cm” s) or a volume or volumetric source strength Sy 
(MeV/cm? s). For example the average gamma source strength Sy (f) in the energy 
group with mean energy E, per unit reactor core volume is defined as 


Sy (B) =k ByN(B), (1.1) 
where 


Py average power density or average power per unit reactorcore volume, 


N (B) energy release in the spectral group with mean energy E, (MeVA), and 

k 3.1.1019 f/ws, a constant which relates the number of fissions in a 

reactor to the energy release. Its numerical value is based on an 
assumed release of 200 MeV per fission. 


In a cylindrical reactor core the thermal neutron flux follows a cosine distri- 
bution along the core centerline and a J,,=distribution along a radius from the center 
of the core. Thus also the reactor power and the gamma ray source strength have 
the same distribution. The relation between the maximum value, at the center of 
the core, of a quantity with this distribution, and its average value over the core, 
can be described as the product of two factors f. and f the radial and axial form 
factors, respectively. Their definitions are 


2.405 
1 1 i 


H,/2 Hy 

ah 2H 


tj tadial coordinate of the i:th fuel element 

h axial coordinate, measured from the center of the reactor core 

n total number of fuel elements (112 in R3b), and R, H and H, as de~ 
fined in table 1. Inserting the numerical values for the R3b, one ob- 


where 


tains 
f, = 1.586 
f, = 1.401 
ff, = 2.22 


In this way the axial source strength distribution of the i:th fuel element, re- 
garded as a line source, can be written 


h 
(h) = Shi (0) cos (1.4 


where 


k fy fg P N(E) 


Spi (9) = n Hy (1.5) 
and P = total power (94 MW for R3b). 
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Inserting the numerical values one obtains 


H cm s 


N fo 


e = 


2. GAMMA ATTENUATION THEORY 


The gamma ray flux from a line source, or in general from any extended 
(line, surface or volume) source, is always computed as a sum or integral over 
the differential point source elements into which the extended source can be sube= 
divided. Thus the basic radiation problem is to compute the flux from a point 
source. When this is known the flux from any extended source is only a geomete- 
rical and mathematical problem, Simple as it may seem, these problems can in 
practice seldom be solved, particularly not by manual methods, Throughout this 
work we must for instance approximate (1.4) by (1.5), a constant, because a line 
source with a strength distribution like (1.4) leads to unsolvable integrals. 

The probability that a photon, which has been emitted from a point source of 
strength S, is found within unit solid angle 142 at a radius r is 


(2.1) 
The probability that the photon will travel a distance r without collision is 
(2.2) 


where 
r 
b = face (2.3) 


and p(r’) is the total macroscopic collision crossesection, also denoted total ab- 


a 


f 
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sorption coefficient u,.;. In the case of n finite homogeneous shields of thickness 
t; and total cross-section w;, eq. (2.3) simplifies to 


b= > « (2.4) 


The uncollided flux ®, at a distance r from the point source, or the so called 
point kernel, is then defined as the product of the source strength S and the prob- 
abilities (2.1) and (2.2) 


5 
(2.5) 


If S denotes the number of photons emitted per second, @,, is the number flux 
(photons/em?s), and OE is the energy flux (MeV/cm?s). In the following we shall 
use a slightly different notation. By S we shall mean the source strength in MeVA, 
so that @, is the energy flux. 

The conditions under which the real flux ® equals the uncollided flux ®,, are 
called narrow beam conditions. The reason for this is evident from a consideration 
of the interaction processes which form the total collision cross-section. 

The most important gamma-ray interaction processes out of about a dozen pro- 
cesses are 
1) the photoelectric effect, 

2) pair production, 
3) the Compton effect. 


1) In the photoelectric effect the photon transfers all its energy to an electron of 
the target atom, which thereby is emitted from the atom. Thus the photon is ab- 
sorbed, and since the electron will not travel far in a shielding material, all the 
gammaeray energy can be said to be absorbed in the immediate vicinity of the 
target atom. 
2) In pair production all the energy of the photon is consumed by the creation of 
an electron and a positron. The same argument holds as in the previous case; thus 
pair production is also an absorptive process, 
3) The Compton effect is essentially a scattering process in which the photon col- 
lides with a free electron and then continues its pathinto a changed direction with 
reduced energy. Thus the photon never disappears, but one can still represent the 
Compton effect by the sum of a scattering cross-section and an absorption cross 
section, if one by absorption means the average fractional energy reduction in the 
scattering collision. 

The total cross«section y+, is for historical reasons called total absorption coef- 
ficient, which is incorrect since it also includes pure scattering. Another cross- 
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section in use is the energy absorption coefficient Men: which equals the total 
cross-section less the pure scattering cross-section. 

Only in a very narrow beam of gamma rays, where even a small angular de- 
viation would bring the photon out of the beam, can a photon disappear entirely 
by a scattering collision and the real flux be represented by (2.5) (where p is the 
total absorption coefficient 

In the general or broad beam case a considerable number of Compton-scattered 
photons undergo one or more scattering collisions in a thick shield and thus have 
a fairly large probability to stay within the beam. This means that the uncollided 
flux (2.5) always represents an underestimate of the real flux. 

In order to find an expression for the total flux @ in analogy with (2.5) one 
would have to consider both angular and energetic distributions of the scattered 
photons in space, but in all geometries save possibly for the simplest (homogene- 
ous one=medium problems) this is a hopeless task, particularly with manual methods, 
One therefore has chosen to represent the total flux by . 


= Bd, , (2.6) 


where B is called “build-up factor", and is a function of photon energy, shielding 
medium, thickness of shield (in dimensionless relaxation lengths b) and sometimes 
angle of incidence on shield, Further B is also a function of the geometry, but 
the customary build-up factors are always computed for infinite or semi-infinite 
media or infinite slabs of finite thickness, so that we have no choice of build-up 
factors for different finite geometries. 

Depending on whether the fluxes in (2.6) are number fluxes, energy fluxes, 
fluxes of absorbed energy, dose rates etc. there are corresponding build-up factors 
in use. In manual calculations it is customary to use dose build-up factors both 
for energy fluxes and for dose rates, because the difference is very small and the 
work thereby may be reduced by a factor of 2. 

Further one can set up relations for other sources than point sources in analogy 
with (2.6). Most common are the build-up factors for isotropic point sources and for 
isotropic or monodirectional plane sources. 

The introduction of the build-up factor B evidently does not solve the problem 
of finding an expression for the total flux ®. We still have to find an expression 
for B, and that is a problem to which there are no general solutions, hardly even 
solutions within a very restricted range of parameters. 

The advantage in expressing the total flux ® in the very simple form (2.6) is, 
that once the B’s are known numerically, (2.6) is very convenient to use. The 
numerical evaluation of the B’s has been going on for a long time now on the 
largest computing machines, using Monte Carlo methods and moments methods 
(cf. Goldstein, 1959), and a thorough compilation on the accumulated data 
can be found in Goldstein’s book, 
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In the same way as the uncollided flux from extended sources is found by in- 
tegrating point kernels (2.5) over the geometrical extent of the source, the total 
flux ®, of extended sources is found by integrating the expression (2.6): 


Dy = dv (2.7) 
v 


over the extent v of the source. Within the range of integration one is therefore 
compelled to find analytical expressions for B. They are always empirical, and 
have to be based on existing numerical data. 

It may sometimes be quite difficult to find useful analytical representations for 
B. In some cases the form 


B=1+b (2. 8) 


represents a good fit to the numerical data, particularly for small values of b. 
When b>1, a better fit would be 
B =b, (2.9) 


which certainly always is incorrect for b<1. 
For manual calculations in one medium a very useful form has been suggested 
by Taylor (1954) 


B = A exp + exp (2.10) 
which contains three parameters A, @, and a independent of b but, of course, 
dependent on source energy E,, material, source geometry, etc. The advantage 


of this form is, that it leaves the point kernel (2.5) unchanged, as can easily be 
shown by direct substitution of (2.5) and (2,10) into (2.7): 


b 


= ab - = 
aC exp (- ab) + (1-A) exp ( 2 dv 
-b’ 
+ (1-4) f (2.11) 


where b’ = (1 + @,)ur, and 
Go) mt. 


This holds evidently only for one medium. The situation is much more compli- 
cated for several media, but this problem lays outside the scope of the present 
work. 


3. CORE APPROXIMATION 1: HOMOGENEOUS, UNIFORM SOURCE 


In the present and the following chapters we want to compute the gamma ray 
dose rate at the cylindrical surface of the R3b reactor core. The simplest approach 
then is to assume the core to be a homogeneous mixture of fuel and moderator of 
density 1.665 g/em®, as given in table 1. 

The axial cosine=distribution of the source strength S, cannot easily be accounted 
for in manual calculations. Since we are interested in the dose rate at a point on 
the “equator” of the cylindrical core, where the axial distribution has its maximum, 
we choose the axial distribution constant, equal to its value in the equatorial plane. 

The simplest approximation for the radial J,,-distribution is also a constant, which 
equals the average value over a radius rather than the maximum value at the center 
of the core. Thus the core constitutes a cylindrical source with uniform strength 


S, =f (3.1) 


from the definitions (2.1) and (2.3). 

Following a suggestion by Glasstone (1955) we simplify the geometry by sub- 
stituting a semi-infinite half-space for the finite cylindrical core. This is permissible 
when the dimensions of the cylinder are large (measured in dimensionless relaxation 
lengths), as shown for instance by Grotenhuis (1959), who has discussed the 
accuracy of the approximation in some detail, 

The flux can be found by integrating (2.6) over the semi-infinite half-space in 
a spherical coordinate system, centered at the detector: When @, is given by 
(2.5) and B is the point isotropic dose build-up factor in Taylor’s form,we obtain 
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2% 
af dg ii sin 040 rar = 
i=1 
Sy Ao 
where 

r radial distance from the detector, 
Ay, Ap = 1-A,, and are build-up constants, 
7 total absorption coefficient for the homogenized core, and 
Sy is given by (3.1). 


Actually the formula (3.2) is an improvement on Glasstones formula for 
three reasons. One of them is evident: Glasstone ) does not use any build-up, 
probably because he considers it an unnecessary complication to the anyway crude 
formula. This is an understatement, however, because the uncollided formula of 
Glasstone, 


(3.3) 


is not conservative and a build-up correction is well needed to increase the anyway 
too small values of (3.3). 

The two other improvements lay in the definition of Sy. Glasstone uses ” 
rather than Sy as defined by (3.1), which is also unconservative. Furthermore, in 
the definition of S,, equation (2.1), Glasstone in an example uses N = Neo 
instead of N = N,,.. This is conservative but erroneous, because one thereby neg- 
lects the shielding effect of the individual fuel element in which the gamma rays 
were born, and one overestimates the flux by a factor of 2.7 in the R3b case. 

The build-up constants for a compound material like the homogenized core 
must be determined empirically from the build-up factor curves (B a continuous 
function of the distance b = wr with the energy as parameter), which in turn can 
be found by using the interpolation methods of Goldstein and Wilkins (1954) 
in the curves of their report. 


The numerical values of the build-up constants used are listed in table 3, 


1) Actually the authors of the shielding chapter in Glasstone’s book are E, P, 
Blizard and H. Goldstein. 


Table 3. 


Taylor build-up constants used in the R3b calculations. 


Energy Homogenized core Pure DgO 

6.45 0.0755 0.0912 3.0 0.0405 0.0916 
4.7 0.0524 0.0245 3.0 0.0647 0.0562 
2.4 0.1169 0.0336 2,2 0.1266 0.0693 
2.0 0.131 0.024 2.0 0.131 0,024 
1.7 0.136 0.028 2.5 0.0911 0.0629 


Inserting the numerical values of N.., (E) from table 2, Py from table 1 and 
k and f into (3.1), the values of S, together with the values of the build-up 
constants of the homogenized core from table 3 give, when inserted into (3.2), 
the expected flux values. The resulting flux values together with the corresponding 
dose rates are listed in table 4. 


Table 4. 


Fluxes and dose rates by approximation 1. 


Energy group Flux Dose rate 
(MeV) (1012 Mev/em?s) (107 rf) 
O- 3 5.48 0.88 
3 5 2.11 0.27 
5 8 0.54 0.06 
O0- 8 (not additive) 1.21 
O- 4 5.42 0.78 
4-8 1.45 0.16 
0-- 8 (not additive) 0.94 


It is evident from the table that the method is sensitive to the choice of energy 
groups. The main reason for this is the large uncertainty in the A, values for low 
energies. If we had chosen E in the energy group 0° 3 MeV slightly different 
from 2 MeV, the resulting dose rate 0.88 - 107 rf would have changed consider- 


ably. 
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4. CORE APPROXIMATION 2: HOMOGENEOUS, NONUNIFORM 
SOURCE 


Approximation 1 can be refined in three basically different ways: 

(1) The semi-infinite half-space geometry could be exchanged to the proper 
cylindrical geometry. This refinement is not expected to improve the situation 
very much because of the large dimensions of the cylinder, but it does complicate 
the calculations very markedly. 

(2) The uniform source distribution could be exchanged to a nonuniform, ap- 
proximately correct radial distribution. This will be called core approximation 2. 

(3) The homogeneous core could be exchanged to the true heterogeneous core 
with distributed sources. This refinement essentially includes both of the previous 
steps, so that it should be introduced after approximation 2. We will call it ap- 
proximation 3. At this point the manual work already becomes so formidable, that 
one has to consider computing machines. 

Approximation 2 consists in replacing the radial J,-distribution by a cosine dis 
tribution or, if the coordinates are reversed (it should be recognized that we still 
have a "slab" geometry and not a cylindrical geometry), a sine distribution, The 
Jgndistribution cannot be used because it gives integrals which contain products of 
Jq (Xq ~ x) and the exponential integral E, (ax + b), which cannot to our knowl- 
edge be solved. 

Even a sine distribution cannot exactly be used, because the integral of 


sin x - E, (ax + b) 


gives an infinite number of terms. This is in effect to replace the sine distribution 
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by its Taylor expansion. Retaining only the first term the flux ® becomes 


i=l 


where R’ =R R,» and R and R, are defined in table 1. 
Sy is here given by 


S, =f, f, Sy. (4.2) 
where f, is not defined by (2.2) but rather by 


2 
=— (cos d +d sind - 1), (4.3) 
d2 
2R 


Substitution of numerical values gives 
f, = 2.26, 
and the results in table 5. 


Table 5. 


Fluxes and dose rates by approximation 2. 


Energy group Flux Dose rate 
(MeV) (101? Mev/cm?s) (10 rf) 
O0-- 3 5.97 0.96 
3 5 2.43 0.32 
5 8 0.98 0.11 
0 8 (not additive) 1.39 
O- 4 6.09 0.87 
8 1.74 0.20 
O-- 8 (not additive) 1,07 


Fs 
- 
Ro 
r co dr = 
f, Ry 2R 
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5. CORE APPROXIMATION 3: DISTRIBUTED LINE SOURCES IN 
A HOMOGENEOUS MEDIUM 


There are two well known methods to transform a cylindrical source S, into an 
equivalent line source 


a) In the self-absorption distance method, the strength of the line source is 
S, = , (5.1) 


where r is the radius of the element, and S, contains the total energy release 
Nrot- The attenuation within the fuel element is accounted for by locating the 
line source a distance z within the fuel element. z can most conveniently be 
found from the graphs of Rockwell (1956). 

In the escape probability method, (5.1) still holds, but S, now contains No... 
rather than Not? so that there is no further need for attenuation within the 
fuel element. Geometrically the line source is located at the axis of the 
cylinder. 


b 


Approximation 3 has been investigated under 3 slightly different circumstances. 
Using the above lettering for the two methods, approximation 3 is subdivided into 
methods 


3ay, 3b; and 3p 


The indices stand for: 

1 each fuel element is immersed into pure DoO. 

2 each fuel element is immersed into a homogeneous mixture of the D,O 
and the fuel elements located between the fuel element considered and the 
detector point. 


In method 3a, a subdivision into 3 energy groups has been used, whereas in 
the methods 3b, and 3b, only two groups have been used. 
The flux in method 3a, is given by 


112s, (0) 


2 
a + > Aj seci{fia + aj) maj + Ben oi}. (5.2) 
=1 j=1 


where 


SLi (O) is given by (1.5) with N(E) equal to N. ote)» 
Zi self-absorption distance of the i:th element, 
aj perpendicular distance from the surface of the i:th element to 


the detector point, 


Aj, a; Taylor build-up parameters for pure at the energy E, table 3, 
for DgO at the energy E 

Hen Men for fuel element at the energy E 

6; = arc tg th 


aj+2j 
seci (b,@) is the modified secant integral 


6 
seci (b, 6) = { exp (-bsec x)dx, (5.3) 


plotted by Rockwell (1956) under the notation f (@,b). 
The flux in methods 3b; and 3b, is given by 


112 


Aj seci [(1 bj, 9; il: (5.4) 


i=1 j=1 


where 


S,; (O) is given by (1.5) with N(B) equal to N... (E), 


b; =a + e, 
1) radius of fuel element, 
Ho 
6; = arc "8 3b,” 


In method 3b, gs; = mw andy, A; and @; are as in method 3a, . 
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Table 6. 
Fluxes and dose rates by the approximations 3. 
Energy group Flux Dose rate 
Method 
(MeV) (101? Mev/em’s) (107 rf) 
3a, 3 3.49 0.56 
3 5 2.73 0.36 
5 8 0.95 0.11 
0+ 8 (not additive) 1.03 
3b, 4 4.53 0.65 
4+ 8 1.81 0.21 
0+ 8 (not additive) 0.86 
3b5 O- 4 3.19 0.45 
4-8 1.15 0.13 
O- 8 (not additive) 0.58 
In method 3b, 
Bj B soe for the homogenized sector enclosed by the core periphery and 


an arc with radius a; centered at the dose point, 
Aj, @; Taylor build-up parameters for the same medium. 
Concerning the numerical values of the parameters, it can be mentioned that 
Zz; = 3.0 cm for all elements for which = > 10, 
seci (b, seci [v. 4] for all the considered angles 6; 


The results of the manual calculations with the expressions (5.2) and (5.4) are 
given in table 6. 
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6. CONCLUSIONS 


In table 7 we collect the dose rate results for all energy groups from the tables 
4, 5 and 6, to facilitate a comparison. Although none of the methods can be as- 
sumed to give an absolutely correct result, upon which a quantitative comparison 
could be based, several qualitative conclusions can be drawn from a relative com- 
parison of the methods, 


Table 7. 


Dose rates (107 r/h). 


eee 2 0.88 0.96 0.56 
3 0.78 0.87 0.65 0.45 
4 0.27 0.32 0.36 
ae 6 0.16 0.20 0.21 0.13 
1 0.06 0.11 0.11 
2 group 0.94 1.07 0.86 0.58 
total 
ae oy 3 group 1.21 1.39 1.03 
total 
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The dependance on the number of energy groups is shown by the total dose rates 
of methods 1 and 2. The 3-group results are larger than the 2-group results by a 
factor of 1.29. In itself this factor does not say that 3-group calculations are con- 
servative compared to 2egroup calculations; it depends. largely on the choice of 
mean energies E, and another choice could well have resulted in a factor 1. 
However, it shows that the division into only two groups is much too coarse, The 
division into three groups is maybe also too coarse, but this cannot be judged from 
the present results. The ultimate criterion is, that a division into n or nel groups 
should not produce significantly different results, and that the width E,,,,/ of the 
groups should be so small, that the choice of mean energies does not influence the 
results significantly. 

If we divide the total (3-group) dose rate of method 3a, by 1.29, we obtain 
0.80 - 107, the probable result of method 3a, if it had been performed for 2 
groups. This figure can be compared with the total dose rate by method 3b,. Although 
the difference is not very significant, it seems that the escape probability method 
(b) is conservative compared to the selfeabsorption distance method (a). Since the 
escape probability approach anyway is more straight-forward and simpler to use, it 
must be recommended as the better method, 

The results of approximation 1 clearly show its limitations. At high energies it 
gives far too low values, because the constant source distribution (3.1) underesti~ 
mates sources far from the detector. At 3 MeV or lower, sources far away are of 
little importance, whereas the distribution (3.1) overestimates sources close to the 
detector. We can conclude that approximation 1 gives conservative values when 
there are only low-energetic sources present, as in the case of a shut down reactor, 
where only the longelived fission product activity is present. 

If we next turn to approximation 3, it is evident from the 3bp data, that the 
homogenization of moderator and fuel elements has a very strong influence on the 
attenuation. At 6 MeV the 3b, value is so unconservative that the method cannot 
safely be recommended. Since the difference between 3b, and 3b, is so significant, 
even at 3 MeV, where most of the dose rate is contributed by the fuel elements 
closest to the detector, it can be concluded that the core medium becomes too 
dense by the method of homogenization used in 3b,. 

Approximation 2, finally, is in general conservative. At high energies, when 
a large number of fuel elements contribute significantly to the dose rate, the co» 
sine source distribution is 4 good approximation to the more realistic methods 3a, 
and 3b,. This is confirmed by the dose rates obtained. At low energies again, 
approximation 2 gives higher dose rates than approximation 3, because the cosine 
distribution assumption introduces sources in an otherwise sourcefree region, viz. 
between the detector and the closest fuel element. 

It can be concluded that 3b, is the best method, if one requires a method 
which contains a minimum of approximations and gives results which are not ex- 
cessively conservative. Furthermore it is well suitable for machine programming, 
because the accuracy can easily be increased by increasing the number of energy 
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groups and by taking account of the neglected axial source strength distribution. 
A program based on 3b, can also easily be generalized to calculate the flux at 
any point within the reactor core and on its surface. 

A program of this kind is Garlic, which has been described in a report by 
the author (1959 b). 
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